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In this paper, we introduce the counterterms that remove the non-logarithmic divergences 
of the action in third order Lovelock gravity for static spacetimes. We do this by defining 
the cosmological constant in such a way that the asymptotic form of the metric have the 
same form in Lovelock and Einstein gravities. Thus, we employ the counterterms of Einstein 
gravity and show that the power law divergences of the action of Lovelock gravity for static 
spacetimes can be removed by suitable choice of coefficients. We find that the dependence 
of these coefficients on the dimension in Lovelock gravity is the same as in Einstein gravity. 
We also introduce the finite energy-momentum tensor and employ these counterterms to 
calculate the finite action and mass of static black hole solutions of third order Lovelock 
gravity. Next, we calculate the thermodynamic quantities and show that the entropy calcu¬ 
lated through the use of Gibbs-Duhem relation is consistent with the obtained entropy by 
Wald’s formula. Furthermore, we find that in contrast to Einstein gravity in which there 
exists no uncharged extreme black hole, third order Lovelock gravity can have these kind 
of black holes. Finally, we investigate the stability of static charged black holes of Love¬ 
lock gravity in canonical ensemble and find that small black holes show a phase transition 
between very small and small black holes, while the large ones are stable. 


I. INTRODUCTION 


An interesting framework for studying the non-perturbative quantum field theories is through 
the use of anti-de Sitter/conformal field theory (AdS/CFT) correspondence l|. According to this 
duality, in principle, one can perform gravity calculations to find information about the field theory 
side or vice versa. In this context, the central charges of the dual theory (CFT) relate to coupling 
constants of its dual gravity. Therefore, Einstein gravity with one coupling constant restricts the 
dual theory to a limited class of CFT with equal central charges jsj]. For extension of the duality 


beyond this limit, one needs to involve higher curvature terms in the gravity action. It is clear 
that each correction term introduces a new coupling constant and therefore one may have CFT 
theory with different central charges. Indeed, this procedure leads to the richness of the CFT 
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theory [3]. The most natural extension of general relativity with higher curvature terms and with 
the assumption of Einstein - that the left hand side of the field equations is the most general 
symmetric conserved tensor containing no more than two-derivatives of the metric - is Lovelock 
theory. Lovelock [4] found the most general symmetric conserved tensor satisfying this property. 
The resultant tensor is nonlinear in the Riemann tensor and differs from the Einstein tensor only 
if the spacetime has more than 4 dimensions. Although Lovelock gravity leads to second-order 
field equations and it has ghost free AdS solution f|, it has been recently shown that quadratic 
and cubic ^gravities entail causality violation and there are stringent conditions on the coupling 
constants 

The problem with the total action of Einstein gravity is that it is divergent when evaluated on 
the solutions |?-9]. Due to this fact, all the other conserved quantities which is calculated through 
the use of this action is also divergent. One way of eliminating these divergences is through the use 


of background subtraction method of Brown and York [7]. In this method, the boundary surface 
is embedded in another (background) spacetime, and one subtracts the action evaluated on the 
embedded surface of the background spacetime from the total action. Such a procedure causes the 
resulting physical quantities to depend on the choice of reference background. Furthermore, it is not 
possible in general to embed the boundary surface into a background spacetime. For asymptotically 
AdS solutions of Einstein gravity, one may remove the non-logarithmic divergences in the action 
by adding a counterterm action which is a functional of the boundary curvature invariants lcj, 


11]. Indeed, this counterterm method furnishes a means for calculating the action and conserved 


quantities intrinsically without reliance on any reference spacetime 


12-14], Although there may 


exist a very large number of possible invariants, only a finite number of them are non-vanishing 
in a given dimension on a boundary at infinity. This method has been applied to many cases 
such as black holes with rotation, NUT charge, various topologies, rotating black strings with zero 
curvature horizons and rotating higher genus black branes [isl . Although the counterterm method 
applies for the case of a specially infinite boundary, it was also employed for the computation of 
the conserved and thermodynamic quantities in the case of a finite boundary 


LJLJLC < 

13 - 


All of the works mentioned in the previous paragraph were limited to Einstein gravity. Although 
the counterterm of Lovelock gravity with flat horizon has been introduced OjOL only a few works 
related to the counterterm method have been done for Lovelock gravity with curved horizon. This 
is due to the fact that even for Einstein gravity, the systematic construction that provides the 
form of the counterterms becomes cumbersome for high enough dimensions [10, 11]. Indeed in this 
method, one should reconstruct the spacetime metric by solving iteratively the field equations in 
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the Fefferman-Graham frame |1 


p 1 

ds 2 = Y^ d P 2 + _ [S{o)ij( x ) + P9{i)ij(x) + P 2 g( 2 )ij{x) + •■•] dx l dx 3 , 


(1) 


4 p 2 ' p 

where g(o)ij(®) is the boundary data of an initial-value problem governed by the equations of 
motion. However, even for Einstein-Hilbert theory, solving the coefficients g( P )(x) in Eq. (JX() as 
covariant functionals of g^o)^) is on ly possible for low enough dimensions. Thus, it is expected 
that the holographic renormalization procedure would be even more complicated in Lovelock grav¬ 
ity because of the nonlinearity of held equations. Indeed, because of the nonlinearity of the held 
equations solving (x) in Eq. ffl as covariant functionals of g( 0 ) (x) would be even more cumber¬ 
some. So, the authors in Ref. [2G] presented an alternative construction of Kounterterms. Instead 
of adding counterterms to cancel the divergence at the boundary explained above, they circum¬ 
vented the difficulties of the standard method by using Kounterterms which depend on the intrinsic 
and the extrinsic curvatures of the boundary. They selected the Kounterterms as the boundary 
terms which are regular on the asymptotic region. Indeed, the regularization process is encoded in 
the boundary terms already presented and there is no need to add further counterterms. 

The exact rotating solutions of Lovelock gravity with curved horizon are not introduced till 
now. Indeed, only static solutions of Lovelock gravity with different matter fields are known 21]. 


So, because of the difficulties of the holographic renormalization procedure in Lovelock gravity 
and the nonexistance of an exact rotating solution of this theory, we limit ourselves to the case 
of counterterms of Lovelock gravity for static solutions. The counterterms of asymptotically AdS 
static solutions of Gauss-Bonnet gravity have been introduced in Ref. 3, 3]- Also, the finite 
action and global charges of asymptotically de Sitter static solutions has been obtained in Ref. 


24]. 


Here we like to apply the counterterm method to the case of the static solutions of the held 
equations of third order Lovelock gravity with curved horizon. We define the cosmological constant 
in such a way that the maximally symmetric AdS spacetime 


ds z = — I k + dt 2 + ( k + 


L 2 J 


L 2 


-l 


dr 2 + r 2 dS^ n _ x 


( 2 ) 


be the vacuum solution of Lovelock gravity. In Eq. ([2]) dT, 2 n _ 1 is the metric of an (n — 1)- 
dimensional maximally symmetric space with curvature constant (n — 1 )(n — 2 )k and volume 
Vk,n-i- Indeed, this choice of cosmological constant makes the asymptotic form of the solutions 
of Lovelock gravity to be exactly the same as that of Einstein gravity. Thus, we expect that 
the counterterm introduced for Einstein gravity in [3] may remove the power law divergences 
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in the action of Lovelock gravity. Although the counterterms which should be added to Gauss- 
Bonnet gravity in order to remove the power law divergences of the action for static solutions are 
introduced in Ref. [22], they depend on the Gauss-Bonnet coefficient. However, because of our 


choice of the cosmological constant, our counterterms are the same as those of Einstein gravity 
and are independent of Lovelock coefficients. In order to check our counterterms, we calculate the 
finite action and the mass of the black hole through the use of counterterm method. Then, we 
use these finite quantities and the Gibbs-Duhem relation to obtain the entropy. We find that the 


calculated entropy of the black holes is consistent with the Wald’s formula 


251 ]. As another test of 


our counterterm method, we show that the mass obtained through the use of counterterm method 
satisfies the first law of thermodynamics. We, also, perform a stability analysis of the black hole 
solutions in canonical ensemble and investigate the effects of third order Lovelock term on the 
stability. 

This paper is organized as follows. In section m we review the well-defined action of Lovelock 
gravity. In section m we introduce the counterterms for third order Lovelock gravity for static 
spacetimes. We also, introduce the finite stress energy tensor of this theory. Section HVl is devoted 
to the thermodynamics of the black hole solutions of the theory. We calculate the finite action, 
the total mass, the temperature, the charge and the electric potential. We calculate the entropy 
through the use of Gibbs-Duhem relation and Wald formula and find that they are consistent. We, 
also, investigate the first law of thermodynamics. In Sec. 0 we investigate the thermal stability 
of the solutions in canonical ensemble. We finish our paper with some concluding remarks. 


II. ACTION AND FIELD EQUATIONS 


a 


The bulk action of Lovelock gravity in n + 1 dimensions may be written as [4] 

[n/2] 


I bulk — , .. [ dx yj ' Otp {dip 2A p ) + I ma t 

iGvr pi 


with [x] denoting the integer part of x, a p ’s (p > 2) are Lovelock coefficients, 


n c^l^2• • *^2p— 1 ^2p T)b\b2 ^2p 

L " P ~ 2P° b i b 2--- b 2 P -i b 2 P n a 1 a 2 '-- n °-2p-ia2p 


( 3 ) 


is the Euler density of a 2p-dimensional manifold, b2 P is the general asymmetric kronecker delta 
and 


A p — 


(—1 ) p n(n — 1 )(n — 2) • • • (to — 2p + 1) 
TL2p ' 


( 4 ) 
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The action ([3]) is written in such a way that the maximally symmetric AdS spacetime ([2]) is the 
vacuum solution of action ©. In this notation, the independent coupling constants are L and all 
the Lovelock coefficients. 

From a geometric point of view the Lagrangian of the action Q in 2[n/2] + 1 and 2[n/2] + 2 
dimensions is the most general Lagrangian that yields second order field equations, as in the case of 
Einstein-Hilbert action which is the most general Lagrangian producing second order field equations 
in three and four dimensions. In the rest of the paper, we work in a unit system with aq = aq = 1 
and the dimensionless Lovelock coefficients a p defined as 


_ _ (n — 2)...(n — 2p + 1) 
a P = £2(p-l) a P ’ 


p > 2. 


(5) 


With the definition ([5]) . the cosmological constant for AdS spacetime is 


, , x [n/2] 

n(n — 1 v-', 

A=-W^E("1^ 


2 L 2 


p =i 


In this paper, we consider the third order Lovelock gravity in the presence of electromagnetic 
field. Thus, the action of matter field is 


Imat = J dx n+1 ^gd [fl A u] d [fl A u] , 

where T, t is the electromagnetic potential. The first term in Lovelock Lagrangian is the Einstein- 
Hilbert term R, the second term is the Gauss-Bonnet Lagrangian £2 = R^sR^^ ~ 4:R flu RP u -\-R 2 , 
the third term is 


r n o p PP T 1 Q p l 11 ' T? crKi PP T 1 0/1 jdH'Vgk p p/9 

Z^3 — ZJrC rCcr^priX ^ + oil ^it ut-^ uk ' Z aix Hcwi/pix t 


ap- LV [LK 




+mR lwrJK R. aNlu + 24 R^ aK R au R KV + 16 R^R va R° n - YIRR^R uv + R 3 


( 6 ) 


and the cosmological constant is 


A = - n{n 2L2 l \ l-a 2 + az). 


(7) 


As in the case of Einstein-Hilbert action, the action Q does not have a well-defined variational 
principle, since one encounters a total derivative that produces a surface integral involving the 
derivatives of Sg /W normal to the boundar y d A4. These normal derivatives of 5g pu can be canceled 
by the variation of the surface action [17|, [26j 


x , [n/2] 

Isur — q / d Xyf—^y % ' Ol p Q p , 
87T JdM pt 
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where 


1 

Qp = pf <**[£:£:$ W * 


x ( - K 2 \ 3 (7) - t*K?K? 

' 2 92.73 v ' 32 33 


- jfop-2i2p-l / X _ t 2 K i2 P -2 K i2 P -l \ /g\ 

2 32p-232p-l v '' 3 2p-2 J2p-1 I V > 


In Eq. dSJ) 7 a fe and A)j& = — 7 a are the induced metric and extrinsic curvature of the boundary 
cbW, respectively. The explicit form of the first three terms of Eq. ([8]) are fill ]: 


4ur = J- f d n Xy/^K 

JdM 


j(2) = 

4tt 

( 3 ) = 803 

sur 8 vr 


'974 


d n x v / -7 ( J - 2(?i 1 6 ) A" afe 


f d n Xy /—7 | J 


d n Xy ^ { P - 2G^K ab + 2RJ 


-12R ab J ab - iRabcd ( 2K ac K°K ea - I<K ac K‘ 


ac tA> rxed 


r CLC T/'bd 


(9) 


where J and A are the traces of 

J ab = ^{2KK ac K c b + K cd K cd K ab - 2K ac I< cd K db - K 2 K ab ) 

and 

P ab = ^{[I< A - 6I< 2 I< cd K cd + 8KK cd K d K ec - 6K cd K de I< ef Kf c + 3 (K cd K cd ) 2 ]K ab 
-( 4 K 3 - 12KK ed K ed + 8K de KjI<^ d )K ac I< b - 2AKK ac K cd K de Kf; 

+( 12 K 2 - 12 K ef K e f)K ac K cd K db + 24 K ac K cd K de K e fK bf }, ( 10 ) 

respectively. In Eq. (|9|) G^ b is the Einstein tensor, R abcd ( 7 ) is the intrinsic curvature and G^ b is 
the Gauss-Bonnet tensor of the metric 'y ab given as 

G { ab = 2 (RacdeR b cde - 2 R adbd R cd - 2R ac R c b + RR ab ) - X -t 2lab . 

III. COUNTERTERM METHOD FOR STATIC SOLUTIONS OF THIRD ORDER 

LOVELOCK GRAVITY 


It is well known that the action I bu ik+Isur is not finite for asymptotically AdS solutions. Inspired 
by AdS/CFT correspondence, one needs to add counterterms to the gravity action in order to get 
a finite action. These counterterms are made from the curvature invariants of the boundary metric 
with the coefficients of the higher curvature terms chosen so that power law divergences in the bulk 
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are canceled for all possible boundary topologies permitted by the equations of motion. At any 
given dimension there are only a finite number of counterterms that do not vanish at infinity. This 
does not depend upon what the gravity theory is - i.e. whether or not it is Einstein, Gauss-Bonnet, 
3rd order Lovelock, etc. Indeed, for asymptotically AdS solutions, the boundary counterterms that 
cancel the divergences in Einstein Gravity may also cancel the divergences in Lovelock gravity if one 
chooses the cosmological constant as in Eq. ((4]) . This is due to the fact that the pth order Lovelock 
Lagrangian C v calculated for the metric © is independent of Lovelock coefficients. That is, the 
different orders of Lovelock action do not mix with each other and one may find the counterterms 
for different orders of Lovelock terms separately. This point makes the calculation easier. Of 
course, the coefficients of the various counterterms for different Lovelock terms will be different, 
depend only on L and will be independent of Lovelock coefficients. Thus, using the counterterms 


of Einstein gravity 


Q, 


we may write the counterterms of third order Lovelock gravity as 

p =i 


(ii) 


where 


7 -ip) __ 
1 ct — 


8ir 


j • 

JdM 




dx n \J —yelp < Ap + B p R + C p ( R a R a b — 




n 


4 (to — 1) 


R 


3n + 2 ~ ~ ab ~ n (n + 2) - 3 ^pabpcdp 

i Up / 1 ^ iiii ti'ab t xn f lV )ii ZJX it -tt a cbd 


4(n — 1) 


16(n — l ) 2 


T777“— 2 —R ah D a D b R - R ab D 2 R ab + 


1 


2 (n - 1 )' 


2(n — 1) 


rd' 2 r + 


( 12 ) 


The coefficients A p , B p , C p and D p in Eq. (fT2l) should depend on L. These coefficients for Einstein 
gravity ( a v = 1 , = 0:3 = 0 ) are |l 0 ] 

(n - 1) L 


Ai = — 


L 


B 1 = - 


2(n — 2) ’ 


C'i = - 


L 3 


2(n — 2) 2 (n — 4) ’ 


D 1 = 


L 5 


(n — 2) 3 (n — 4)(n — 6 ) 
(13) 


One may note that the coefficients A pi B p , C p and D p are independent of Lovelock coefficients. We 
apply the counterterms (fl2l) to various static solutions of Gauss-Bonnet and third order Lovelock 
gravity with different topolgy and find that these coefficients are 


A 2 = --^, 


B 2 = Bi, C 2 = - 6 C 1 , 


Z?2 = — 10Di, 


(14) 


and 


^3 = 7^1, 
5 


£3 = ~Bi, C 3 = -18(71, 


£3 = 15£i 


(15) 
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for Gauss-Bonnet and third order Lovelock gravity, respectively. The reason that we use exactly the 
counterterms of Einstein gravity is as follows. First, the boundary at r = const, for static solutions 
is a constant curvature hypersurface and therefore no six-derivative term will be appeared in the 
counterterms. In other words, all the terms of a specific order of counterterms [for example R 2 
and R ab R ab in C 2 (R ab R a b ~ nR 2 /[4(n — 1)])] for static solutions are proportional to r~ 4 . Also, 
Rabcd,R abcd is proportional to r^ 4 . Therefore, in order to remove the divergences of the action 
which are proportional to r -4 , any combination of R 2 , R ab R ab and R abc dR abcd can be used. So, 
we just use exactly the counterterms of Einstein gravity. Second, A p , B p , C p , D p , for p = 2 and 
3 are proportional to those of Einstein counterterm independent of the dimensions. That is, the 
dimensional-dependence of these coefficients are the same as those in Einstein gravity. Third, as 
we will see in the next section, the entropy calculated through the use of Gibbs-Duhem relation 
and the mass and action calculated by our counterterms is consistent with the entropy obtained 
by use of Wald’s formula. Fourth, the mass calculated by our counterterms satisfies the first law 
of thermodynamics. 

While the total action hulk + hur + Rt is appropriate in grand-canonical ensemble where 5A p 
is zero at the boundary, the appropriate action in the canonical ensemble where the electric charge 
is fixed is [27] 

I = hulk + hur + ht ~ tx [ d n x y /^n ,1 F pi/ A u . (16) 

16-7T J dM 

Thus both in canonical and grand-canonical ensemble, the variation of total action about the 
solutions of the field equations is 

3 ( x t(p ) x t(p ) \ 

w = = + < 17 > 

So, the energy-momentum tensor can be written as: 

9 3 / A A \ 

rp _ rr(sur) , rr(ct) _ 4 I 01 sur 01 ct \ 

Tab ~ T ab + T ab - — 2^ J • (18) 

The explicit expressions of and are somewhat cumbersome so we give them in the 

Appendix. 

To compute the conserved charges of the spacetime, we choose a spacelike hypersurface B in 
dAi with metric aij, and write the boundary metric in ADM form: 

= — N 2 dt 2 + dij {dip 1 + V l dt ) (dpi + V^dt) , 


7 ab dx a dx a 


( 19 ) 
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where the coordinates p 1 are the angular variables parameterizing the hypersurface of constant r 
around the origin, and N and V 1 are the lapse and shift functions, respectively. When there is a 
Killing vector field ? on the boundary, then the quasilocal conserved quantities associated with the 
stress tensors of Eq. (fl8l) can be written as 

Q(<t) = [ d n ~ x p\foT ab n a ^ b , (20) 

JB 

where a is the determinant of the metric a t j, and n a is the timelike unit normal vector to the 
boundary B. In the context of counterterm method, the limit in which the boundary B becomes 
infinite (Boo) is taken, and the counterterm prescription ensures that the action and conserved 
charges are finite. No embedding of the surface B into a reference of spacetime is required and the 
quantities which are computed are intrinsic to the spacetimes. 


IV. THERMODYNAMICS OF ADS CHARGED BLACK HOLES 


The field equation of third order Lovelock gravity for the static metric 

ds 2 = — f (r)dt 2 + + r 2 dY ? Kn _ t 

in the presence of electromagnetic field may be written as 

mJ 2 9 n 2 T 2 

*> - 1] - - 1] + «r) - 1] + — - ( „_ 1)( „_ 2)r 2 („-i) = °> ( 21 > 

where ip(r) = [f(r) — k\L 2 /r 2 , q is the charge of the black hole and m is the integration constant 
which is related to the mass of the solution. The electromagnetic potential for the above metric is 


\ —_ 1 & 

M (n — 2)r n_2 ° M ' 

The mass parameter m may be written in terms of horizon radius and q by using the fact that 
i/?(r- 1 _) = —kL 2 /ro_, and therefore 

2 q 2 


L 2 \ rp 

m = I 1 - d 2 + d 3 + k-^fi ) -T + 


r\ ) L 2 ( n — l)(n - 2)? 


n—2 > 


where 


L 2 L 4 

/r = 1 + dt^k— p—b d^k 2 — J-. 


( 22 ) 
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As one expects, ij>(r ) = 1 ( f(r) = 1 + kr 2 / L 2 ) is the root of Eq. (1211) for the vacuum spacetime 
(m = 0 and q = 0 ). The three solutions of the cubic equation (12TT) are 


where 


5 = (v + \J v 2 — u 3 ) 
a 2 — 3d3 


ipi(r) = + <5 + u5 \ 

^(r) = ^ + + 


1/3 


(23) 

(24) 

(25) 


u = 


9 «i 


1 2d 3 — 9d 2 d 3 — 27d2ci3 1 (mL 2 

2 54 d 3 2 d 3 \ r n 


2L 2 q 2 


(n — l)(n — 2 )r 2 ( n_1 ) 

All of the above three roots could be real in the appropriate range of a 2 and 03 . The second and 
third solutions are real provided it 3 > v 2 . Here, we will consider only the first solution ipi(r) which 
is real provided u 3 < v 2 at any r. 

Now, we investigate the thermodynamics of the black hole solutions. The temperature of the 
event horizon may be calculated through the use of analytic continuation of the metric. One obtains 


where 


T = idb ( [1 " “ 2 + “ 3]n & + Kn » ~ 2Tl) ~ (re — l^ 2 ( n ~ 2 ) 


kL 2 k 2 L 4 

f? = 1 + 2d 2 —+ 3d 3 — 3 —. 

7 + r\ 


(26) 


(27) 


The charge of the black holes per unit volume can be calculated by integrating the flux of the 
electric field as 


1 f* 1 

Q = — / FdQ = —q. 

4-7T J 47T 

The electric potential <b, measured at infinity with respect to the horizon is defined as 


(28) 


d* — ^4/iX^lr^-oo ^4/iX^|r=r+- 


(29) 


where x = d/dt is the null generator of the horizon. One obtains 

q 


<f> = 


(n — 2 )r 


n —2 ' 


(30) 
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The finite total action in grand-canonical and canonical ensembles can be found through the 
use of the counterterm method introduced in the last section. It is a matter of straightforward 
calculations to show that the total action is finite. Since we are interested in the stability of the 
solutions in canonical ensemble, we calculate the finite Euclidian action per unit volume 14,n-i i n 
this ensemble. One obtains 


/ = 


16n 




+ 


1 


(n — 1)77 (n — 2 )) r 


^2 + ( n “ 1)^+ 2 


-(1 - d 2 + < 53 ) 


n ( - — 1 ) + 1 
V 


j2 ~ ~ 2 ri)- r + 2 

Tj 


+ h 


where /3 is the Euclidean time period (the inverse of temperature), is given in Eq. (1271) . £ is 

(n — l)kL 2 (n — l)fc 2 L 4 

£ = 1 + 2<5 2 -^--77— 9 —I - 3a s 


and In is 


/n = 


P_ 

8tt 


E (-T 


(n — 3 )r 2 

l/2) ((n - l)!!) 2 L n ~ 2 
n! 


4 > 


(n — 5)? v 


, 0 (« - !) - , 0 (™ - !) - 
1 + 2^- — ot 2 + 3 t -^0:3 


(31) 


(n-3) 


(n — 5) 


n/2=3,4,.. 

One should note that Jo appears only in odd dimensions (even n). 

For our static solution, there is a Killing vector field <, = d/dt on the boundary, and therefore 
the quasilocal conserved quantity associated with the stress tensors of Eq. (fTHl) is the mass of 
the black hole. Using the counterterm method introduced in the last section, it is a matter of 
calculation to obtain the mass of black hole as 


, , in — 1 )m 
M = v ; -1- M 0 



,L 2 

1 — d 2 + d 3 + k —n 


+ 


I 67 t 2 Q 2 


+ M 0 , 


(32) 


L 2 ' ( n - l)(n-2)rl~ 2 ) 
where Mq = Iq/P is the Casimir energy for the vacuum AdS metric per volume 14 n _i, which is 
nonzero only in odd dimensions. In order to have positive energy M — Mo, one should restrict the 
range of Lovelock coefficients in terms of Q , L and r+ as 


0:3 > — ( 1 + k 


Z 6 


-1 


l 2 

^2 


1 + k —^—b <5 2 ( 1 — k 2 


2 r 2 


2 q 2 L 


(33) 


(n — l)(n — 2)r^ n 1 -* J 

That is, the third order Lovelock coefficient has a lower limit in terms of <5 2 , q, L and r + . 

Using the mass and action per unit volume 14 jn _i calculated through the use of counterterm 
method, one may calculate the entropy per unit volume through the use of Gibbs-Duhem relation 
S = pM- I as 


n— 1 


S = £- 


4 


(34) 
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It is worth to note that the entropy of the black hole solution per unit volume T4,n-i calculated 
by the use of G: 
formula fl 28] 


by the use of Gibbs-Duhem relation is consistent with the calculation through the use of Wald 

nn 


S = ^jd n + kR + k 2 a2C 2 ), (35) 

where the integration is done on the (n — l)-dimensional spacelike hypersurface of Killing horizon, 
g^ w is the induced metric on it, g is the determinant of g and £2 is the 2nd order Lovelock 
Lagrangian of g^ u . Also, one may note that the thermodynamic quantities calculated in this 
section satisfy the first law of thermodynamics 


<1M 


(dM/dr+) Q 

(dS/dr+) 


dS + 



TdS + $dQ. 


V. STABILITY IN THE CANONICAL ENSEMBLE 


Now, we study the thermal stability of the black hole solutions of third order Lovelock gravity in 
canonical ensemble. First, we investigate the conditions of having black hole solution in Lovelock 
gravity. The solution given by Eq. (|23j) presents a black hole solution provided f(r) has at least 
one real positive root. This occurs if r + > r ext , or q < q ext where r ext and q ex t satisfy the following 
equation: 

[l-& 2 + d 3 ]n^f + k[(n - 2) + fc(n - 4)^ + A: 2 (n - 6)^-( 36 ) 

^ ext Text (n — 1 )v ex £ 

Also, the temperature of a physical black hole should be positive. One may note that the tem¬ 
perature changes its sign at the root of rj = 0 (r crit ). The radius r crit depends on L and Lovelock 
coefficients, while r ex t depends on q too. In Fig. |H the vertical line is g = 0 line. This figure 
shows that one may have only small (77 < 0 ) and large black holes, and there is no medium black 
hole solution ( r cr it < r + < r ex t). This feature does not occur for black holes of Einstein gravity 
or Lovelock gravity with positive Lovelock coefficients, since r/ 7 ^ 0. Thus, for the case of Lovelock 
gravity with negative Lovelock coefficient(s) one may divide the black holes into two classes with 
negative and positive rj. 

In the canonical ensemble Q is fixed and therefore the black hole solutions are stable provided 
(d 2 M/dS 2 )Q > 0 in the range that T is positive. Using the expressions for mass and entropy given 
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in Eq. 


and (|3ift . one may calculate ( d 2 M/dS 2 )Q: 

fd 2 M\ f dS 
q Vdr 2 ' ' 


d 2 M\ 
'dS 2 ) ( 


Q 


dr .i 


-2 


dM\ 

) Q 


dr 


dS 
dr. l. 


-3 


d 2 S 

dri 


(n — 1 )a + 2 q 2 qr + 
(n — 1 ) 2 r” + 4 ry 3 


—2n+8 


(37) 


where 


L 2 L 4 

7 = ( 2 n — 3) + 2kct2{2n — 5)— + 3k 2 ds(2n — 7) —, 

r + r + 

a = nr®(l — 07 + d^)L ~ 2 + fcr 4 [2 + n (—1 — 60:2 + 6d2(l + da))] 

—fc 2 r^ [n(d 2 ( 15 d !3 + 1) — 15 d 3 (l + 0:3)) — 80:2] L 2 — 2k [n(d 2 — 253) — 4d2 — 60:3] L 4 
~3k 2 dsct2r^ 2 (n — 8 )L 6 — 3/cd|?’^ 4 (?z — 6)L 8 . 



FIG. 1: 10 2 T versus r+ for q = 0.1, L = 1, k = 1 and n = 6 in Einstein gravity (solid) and Lovelock 
gravity (dotted) with 02 = 0.1 and 03 = —0.01. The vertical line is 7 = 0 line. 

To investigate the stability of black holes of third order Lovelock gravity in canonical ensemble, 
we should find the sign of (0 2 Af/55 2 )g, when T is positive. In order to investigate the stability, 
we plot T and (d 2 M/dS 2 )g in one figure. The allowed region for investigation of (5 2 M/95 2 )g 
is when T is positive. In Fig. [2j the vertical dotted-line is 77 = 0 line. This figure shows that the 
small black holes (77 < 0) divided into stable and unstable black holes. That is a Hawking-Page 
phase transition exists for small black holes between very small and small black holes. There is 
no phase transition for large black holes (77 > 0) as one may see in Fig. [2j For investigating the 
effect of third order term of Lovelock gravity on the stability of the black hole solutions, we plot T 
and (3 2 M/95 2 )g versus 0:3 for small and large black holes. Figure [3] shows that there is no large 
black holes for large negative 0 : 3 , while the large black holes for small negative 013 are not stable. 
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FIG. 2: 10 3 T (dotted) and 10 5 (< d 2 M/dS 2 )Q (solid) versus r + for a 2 = 0.1, = —0.01, <7 = 0.1, L = 

1, k = 1 and n = 6 . The vertical line is 77 = 0 line. 



FIG. 3: 10 2 T (dotted) and 10 1 (d 2 M / d S 2 ) ^ (solid) versus 0:3 for a 2 = 0.1, q = 0.5, L = 1, fc = 1, n = 6 
and r+ = 1.2. The vertical line is rj = 0 line. 

As 0(3 becomes larger, large stable black holes may exist. Figure 0] shows that there is no small 
black hole solution for positive 0:3, while the small black hole with negative 0:3 are stable. Figure 0 
shows that there is no negative-?] (small) uncharged black hole with event horizon, while the large 
black holes (77 > 0 ) show a phase transition between large and very large ones. One may see in 
Fig. [5] that the large black holes with positive rj show a phase transition, while these black holes 
with negative rj are not stable. Finally, Figs. [6] and [7] show the temperature and (d 2 M/dS 2 ) q 
versus 0:3 for small and large uncharged black holes. Figure [6] shows that there is no uncharged 
black hole with negative 77, while the black holes with positive r/ are stable. On the other hand, 
the large uncharged black holes with positive rj show a phase transition as <53 becomes larger and 
these solutions are not stable for very large negative 0:3. 
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FIG. 4: 10 4 T (dotted) and 10 4 (d 2 M/dS 2 ) Q (solid) versus a 3 for a 2 = 0.1, q = 0.5, L = 1, k = 1, n = 6 
and r+ = 0 . 1 . 



FIG. 5: 10 2 T (dotted) and ( d 2 M/dS 2 )g (solid) versus r+ for d 2 = 0.1, 0:3 = —0.2, q = 0, L = l, k = l 
and n = 6 . The vertical line is 77 = 0 line. 

VI. CLOSING REMARKS 


The concepts of action plays a central role in gravitation theories, but the sum of the bulk 
and the surface terms diverges. In this paper, we introduced the counterterms that remove the 
non-logarithmic divergences of static solutions of third order Lovelock gravity. We did this by 
defining the cosmological constant in such a way that the AdS metric (|2|) is the vacuum solution of 
Lovelock gravity. Indeed, the cosmological constant J4|) makes the asymptotic form of the solutions 
of Lovelock gravity to be exactly the same as that of Einstein gravity. Thus, we employed the 
counterterms introduced for Einstein gravity in Q and found that the power law divergences of 
static solutions in the action of Lovelock gravity can be removed by suitable choice of coefficients. 
We found that the counterterms are independent of Lovelock coefficients and the dimensionally 
dependent of them is the same as those of Einstein gravity. The main difference of our work with 
that of Ref. [16] is that the counterterms are exactly the same as those of Einstein gravity and 
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FIG. 6 : 10T (dotted) and 10 (d 2 M/dS 2 ) (solid) versus 03 for a 2 = 0.1, q = 0.5, L = 1, k = 1, n = 6 and 
r+ = 0.2. 



FIG. 7: 10 2 T (dotted) and 10 2 (d 2 M/dS 2 ) Q (solid) versus 0:3 for oh = 0.1, q = 0.5, L = 1, k = 1, n = 6 
and r + = 1.2. The vertical line is 77 = 0 line. 

do not depend on Lovelock coefficients. The only job which remains is that one needs to calculate 
the coefficients A p , B p , C p and D p . For example if one wants to remove the divergences of the 
action due to the Gauss-Bonnet term, one should calculate A 2 , B 2l C 2 and D 2 that remove the 
divergences of y/^L 2 , without regarding the other Lovelock terms. This enables one to generalize 
this method to other higher curvature theories of gravity easily, including fourth order Lovelock 
gravity or f(R) gravity. We also introduced the finite energy-momentum tensor in third order 
Lovelock gravity. 

In addition, we employed these counterterms to calculate the finite action and mass of the static 
black hole solutions of third order Lovelock gravity. Calculating the temperature, the electric charge 
and electric potential and using the calculated finite action and mass, we showed that the entropy 
calculated through the use of Gibbs-Duhem relation is consistent with the calculated entropy by 
Wald’s formula. We, also, showed that the conserved and thermodynamic quantities satisfy the first 
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law of thermodynamics. Finally, we investigated the stability of charged black holes of Lovelock 
gravity in canonical ensemble. We found that the black holes with respect to the sign of i] given 
in Eq. ()27l) may be divided into two classes. The negative i] (small) black holes show a phase 
transition between very small and small black holes, while the large black holes are stable. There 
is no black hole solution with medium size. Of course by small black holes we mean the size of 
them with respect to the cosmological parameter L. We, also, investigated the effects of third order 
Lovelock term on the stability of the solutions. We found that there is no large black holes for 
large negative 013, while the black holes for small negative 0:3 are not stable. We also found that 
as a 3 becomes larger, stable black holes may exist. This shows that there is a phase transition as 
the third order term of Lovelock gravity becomes larger. Finally, we found that there is no small 
black hole solution for positive a. 3 , while the black hole with negative a 3 are stable. Finally, we 
considered the uncharged black holes of Lovelock gravity and found that negative t] solutions are 
not black holes with event horizon, while the positive 77 black holes show a phase transition. 
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VII. APPENDIX: 

The first term of Eq. (|18l) which is the energy-momentum tensor of the surface term is 

T ab Ur) = ^{ R ab - K lab + 2a 2 [3J a b ~ Jlab ~ 2 G {a c K b)c + 2 R ab K - K ab R + 2K cd G cd lab - 2K cd R acbd \ 
+3a 3 [5P ab - P lab + 2KG® + C 2 (K ab + K lab ) + 4 JR ab - 24 J {a c R b)c + 8 K cd R ac R bd - 8 KK a c K b d R cd 
+8 KK ab K cd R cd - 8 K cd R ab R cd + l 6 K {a c R b) d R cd + 16 K {a c K b) d K d e R ce - 8 K ab K c e K cd R de + 6 J ab R 
-8 K {a c R b)c R - 12 J cd R acbd - 4 K cd RR acbd - 8 K a c K bc K de R de + 16 K c d K e f K {a c R b)edf + 16 KjK de K\ a c R b)ecf 
+16 K {a c R b)dce R de - 16 KK de K {a c R b)dce - 16 K cd R {a e R b)cde - 4 K cd R ac e f R bdef + 8 K cd R c e R aebd 
+8 K cd R c e R adbe - 8 K cd R a e jR bedf - 8 K\ a c R b) de f R cdef + 8 K a c K b d K e f R cedf - 4 K ab K cd K e f R cedf 
+8 K cd R a e b fR ced} + 2 lab (G^K cd + 6 J cd R cd - JR + 2KK cd I< e f R cedf - 4 K c e K cd I<f h R dfeh )]} 
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and the second term corresponding to the counterterm is 
3 

Tab ^ = g^: °P | A>7afe — 2 Bp (flab ~ ^7-^ 

P= 1 ^ ' 


-a 


- lab ( R cd R cd - 


n 


4 (n — 1) 


R 1 2 I - 


n 


(n-1) 


RR-ab 


1 7afe R 2 R + R 2 R a6 - 44 Da D b R + AR acbd R cd 


—D 


(n-1) 

"2(3n + 2) 
4(n — 1) 


(n-1) 

G^R cd R cd ) - D a D, 


+labD c D d ( RR cd )+D z ( RR ab ) - D c D b ( RR ac ) - D c D a ( RR bc 


b (R e/ R e/ ) + 7 abD 2 (R ef R ef ) + 2 RRZR bc 

n(n + 2) 


-7 abR 3 


+6 ( R~R ab — D a D b R + 'y ab D^R' 


- 2 


16(n — l) 2 

- lab R ei R cd R ecfd + 3 ( R a e R cd R ecbd + R b e R cd R ecad 


-2 D c D d (R a& R cd ) + 2 D c D d ( R a c R b d ) + 2 lab D e D^ ( R cd R e cfd ) + 2D 2 (R cd R acM ) 

—D a RD b R + R b c D c D a R 


-2D e D n ( R cd R e M ) - 2 D,D h ( R ca R 


' cbd 


hcd, fye 


' cad 


+ 


n — 2 
2 (n - 1) 


+R a c D c D b (R) - 2 D c D c D b D a (R) - - R a c b d D c Dd{R) + 2 D a {R bc )D c (R) + D b (R ac )D c {R) 
+D b R ac D c R - 2 D c R ab D c R - R b cd D a D d R - R a cd D d D b R 
- R a c b d D d D c (R ) + 2 R ab D e D e {R) + ^ 7ab (^D e RD e R + D A R 

- 2 D a D b {D 2 R) - 2D 4 R ab + 2D p R (a R P b) - 4 R pq D a D b R pq + 12 R pq D (h R pa)q - 4 {D 2 R pq )R paqb 
+6 D p (R)D {b R P a) - 2D a (R pq )D b (R pq ) + AD {a D q (R pq R b)p ) + 16£> r ( R^I^R^a) 

-4 R pq R pa R qb + 8 R pr R q R paqb + 4 R^R[ a R| r9p , ft) 


+ 7a& 


- R 4 R - 4R p R g (R)RW + 2R pq D 2 (R pq ) + 4 D r D s {R pq )R prqs - D e (R)D e {R) 


+5 D r (R pq )D r (R pq ) -8 D r (R pq )D q (R pr ) - A R pq R p R qr + A R pq R rs R prqs 


+ 


1 


n — 1 


2 D a D b (D 2 R) - 2 (D 2 R)R ab + D a (R)D b (R) + lab ( -2 R 4 R - J D C RD C R 


where Tip’s, Rp’s, Cp’s and D p ’s are given in Eqs. (1131115|) . The above counterterms remove the 
divergences of the energy-momentum tensor for n < 8. As in the case of Einstein gravity, one 
should add more counterterms for n > 8. 
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